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Abstract 

The radial structure of the continuous spectrum of shear Alfven waves is calculated in this 
paper within the separatrix of a magnetic island. Geometrical effects due to the noncircularity 
of the flux surface's cross section are retained to all orders. On the other hand, we keep 
only curvature effects responsible for the beta-induced gap in the low-frequency part of the 
continuous spectrum. Modes with different helicity from that of the magnetic island are 
considered. The main result is that, inside a magnetic island, there is a continuous spectrum 
very similar to that of tokamak plasmas, where a generalized safety factor q can be defined and 
where a wide frequency gap is formed, analogous to the ellipticity induced Alfven eigenmode 
gap in tokamaks. The presence of this gap is due to the strong eccentricity of the island cross 
section. The importance of the existence of such a gap is recognized in potentially hosting 
magnetic- island induced Alfven eigenmodes (MiAE). Due to the frequency dependence of the 
shear Alfven wave continuum on the magnetic-island size, the possibility of utilizing MiAE 
frequency scalings as a novel magnetic-island diagnostic is also discussed. 

PACS: 52.55.Tn, 52.35.Bj 

1 Introduction 

Shear Alfven waves (SAW) are electromagnetic plasma waves propagating as transverse waves 
along the magnetic field, with the characteristic Alfven velocity va = B/^/4ng (B is the magnetic 
field and g the mass density of the plasma). In fusion plasmas, fast ions in the MeV energy range 
have velocities comparable with the typical Alfven speed. In addition, SAW group velocity is 
directed along the magnetic field line and, therefore, fast ions can stay in resonance and effectively 
exchange energy with the wave [TJ [2] . SAW in a nonuniform equilibrium experience collisionless 
dissipation (continuum damping 0111 E])> due to singular structures that are formed where the 
SAW continuum is resonantly excited. In tokamaks, the magnetic field intensity varies along the 
field line. This creates gaps in the SAW continuous spectrum [6] due to translational symmetry 
breaking, analogous to electrons traveling in a periodic lattice [HE]- 

Two types of collective shear Alfven instabilities exist in tokamak plasmas: energetic-particle 
continuum modes (EPM) [7\, with frequency determined by fast particle characteristic motions, 
and discrete Alfven eigenmodes (AE), with a frequency inside SAW continuum gaps [8j. The 
former can become unstable provided the drive exceeds a threshold determined by the continuum 
damping absorption; the latter, have a generally lower instability threshold, being practically 
unaffected by continuum damping [TJ [2J [3J [H [5] . For this reason, the importance of understanding 
the continuum radial structure is clear, if one faces the problem of reaching the ignition condition 
with magnetically confined fusion plasmas, due to the potential impact of AE stability on the 
plasma confinement. 
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The nonlinear dynamics of SAW is a research topic still rich of open issues. In the special case 
of uniform ideal incompressible plasmas a peculiar state exists, called the Alfvenic state [U| I10L fTTj . 
where the nonlinear effects due to Maxwell and Reynolds stresses cancel and give a self-consistent 
nonlinear state. On the other hand, in tokamak plasmas, nonuniformities and non-ideal effects 
such as resistivity or finite plasma compressibility are present. This breaks the Alfvenic state 
conditions and therefore SAW are characterized by a rich nonlinear dynamics. As a consequence, 
the SAW continuous spectrum can be modified by the interaction with low-frequency MHD 
fluctuations [12], such as magnetic islands. Magnetic islands are the result of magnetic-field 
reconnection processes, which occur in tokamak plasmas because of non ideal effects, such as 
finite resistivity [13] . The typical island oscillation frequency and growth rate are much lower 
than the SAW frequency. Due to this time scale separation between island and the SAW dynamics, 
the island is modelled here as a static helical distortion of the equilibrium. 

We derive our fluid theoretical description of the SAW continuum structure inside a finite size 
magnetic island in finite-/? tokamak plasmas. The magnetic island is modelled here as a straight 
flux tube with a noncircular cross section. Therefore, we neglect toroidicity effects on coupling 
modes with different poloidal mode numbers. However, we do take account the low-frequency gap 
in the SAW continuum due to the effects of geodesic curvature and finite compressibility |14[ fT5j 
16J. In our model, the coupling of SAW and the island magnetic perturbation is a linear problem 
of SAW dynamics in a modified static equilibrium. The local differential equation yielding the 
SAW continuum structure is solved numerically with a shooting method code, in the spatial range 
of interest inside the magnetic island [T7] . We are interested here in modes with different helicity 
from that of the magnetic island; in other words, we assume d/dC, 7^ 0, where £ is the coordinate 
along the flux tube. The solution of this problem inside and outside the magnetic island, for 
modes with the same helicity of the magnetic island, is given in a different paper [18) , 

As suggested by the magnetic field line helicity behavior [19} 112] . in an equilibrium with 
a magnetic island a generalized safety factor q can be defined for each flux surface, as an ap- 
propriate average of the ratio of the magnetic field components along and perpendicular to the 
magnetic axis. This safety factor is a monotone function of the distance from the O point, and 
grows towards the separatrix, where it has a singularity. This is due to the vanishing of the 
magnetic field intensity near the X point of the magnetic island. Within the separatrix, we find 
a continuous spectrum very similar to that of tokamak plasmas, where the coupling of modes 
with different mode numbers creates frequency gaps. In particular, a wide frequency gap is 
formed [T7], analogous to the ellipticity induced Alfven Eigenmode [20] (EAE) gap in tokamaks. 
The presence of this gap is due to the strong eccentricity of the island cross section. In this gap, 
a magnetic-island induced Alfven eigenmode (MiAE) could exist as a bound state, essentially 
free of continuum damping, and interact nonlinearly with the magnetic island, if the thermal 
or energetic component of the plasma provide sufficient free energy for driving the mode. Due 
to the frequency dependence of the shear Alfven wave continuum on the magnetic-island size, 
the possibility of utilizing MiAE frequency scalings as a novel magnetic island diagnostic is also 
discussed. 

The scheme of the paper is the following. In Sec. [2] we describe the equilibrium magnetic 
field, starting with the tokamak coordinates and introducing the cylinderlike coordinates for the 
magnetic island flux tube. The generalized safety factor defined inside the magnetic island, is 
shown to have similar behavior of the safety factor in diverted tokamaks. In Sec. El the SAW 
model equations are introduced and the equation for general shear Alfven modes is provided. The 
equation for continuum modes is then derived in the cylinderlike coordinate set defined inside the 
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island and the eigenvalue problem is established. The resulting continuous spectrum, obtained 
by solving numerically the eigenvalue problem with a shooting method code, is shown to be 
very similar to the continuous spectrum of tokamaks. The result is provided for typical tokamak 
parameters, and typical size magnetic islands. The analogy with the tokamak SAW continuum 
is shown in the limit of very large magnetic islands, when the section of the magnetic island 
flux tube is nearly circular. Finally Sec. U] is devoted to a discussion about the applications of 
our results, in gaining a better knowledge on the stability of tokamak plasmas. The differential 
operators calculated in the cylinderlike coordinate set are provided in Appendix |A"1 



2 Equilibrium and coordinate system 
2.1 Tokamak coordinates 

We consider a tokamak geometry where Rq is the major radius of the torus. The equilibrium 
is made of an axisymmetric tokamak magnetic field, plus a helical perturbation generating the 
magnetic island. We start our analysis by describing the axisymmetric field with an orthogonal 
coordinate system (tt,0t,Ct) where tt is the radial coordinate, 9t is the poloidal angle, corre- 
sponding to the field B po i = B po i(rT) Ot, and Ct the toroidal angle, corresponding to the toroidal 
field Btor = Bt or C,Ti considered uniform for simplicity. The subscript T denotes here the tokamak 
coordinates. The gradients associated with these coordinates are: Vry = tt, V#t = 6t/tt and 
VCt = Ct/Ro- 

We want to describe with a slab model, the region of the tokamak plasma where a magnetic 
island with poloidal and toroidal numbers (mi S i,ni S i) is located. This is the region around the 
magnetic island's rational flux surface qT = Qo = fTHsl/ n iaU corresponding to ry = ro, where 
the axisymmetric field has value Bq = jBt or , with 7 = \fX-\- £q/Qq- Here qx is the tokamak 
safety factor qx = rTB tor /(RoB po i), and So = tq/Rq is the inverse aspect ratio, assumed to 
be small: £0 ^ 1- Therefore, we introduce an orthogonal coordinate system (gy,it, £) through 
a rotation in the (#t,Ct) plane, defining u = n^i^T — qo^r) and £ = (Ct + £q &t / qo) / (qol 2 ) ■ 
This is appropriate to describe a sheared field problem, in the proximity of the rational surface, 
because the axisymmetric magnetic field is directed only along £ at qx — qo = 0. The gradients 
associated with these coordinates are: V(/t = (<?os/ r o)^T> ^ u = n isi(CT — qo&T / £0) / Ro and 
V£ = (Ct + £o &t / qo) / {qoRol 2 ) ■ The parameter s is the magnetic shear calculated at the rational 
surface. It is useful to rewrite the gradient of the coordinates u and £ in the form Vu = u/po 
and VC = where po = ro/((/o n isZ7)) and Zq = jqoRo. In fact, we see that 27rpo is the 

length of the magnetic island in the u direction and 2ttZq is the length of a magnetic field line of 
the axisymmetric equilibrium at qT = qo- In other words, 2ttZq is the periodicity length of the 
magnetic flux tube defined as the region inside the magnetic-island separatrix. 

The axisymmetric magnetic field components can be written now by considering that £ forms 
an angle ao = ar ct an (B po i(qo)/ B tor ) with C,t- Therefore in these variables the magnetic field of 
an equilibrium with a magnetic island is B = B qTtP hfT + B UtP hU + B^ tP hC, with: 

d d • 1 \ d £oBo{qr - qo) D D 

Hq T ,ph — £>isl Sm{U), D u ph — 5— , t!(- p h — Dq 

%T 

where we dub the components rescaled with the length of the corresponding basis vector as 
physical components (see Appendix [XJ. We have adopted the constant-"^) approximation for the 
magnetic island, assuming a constant magnetic field amplitude of magnetic island perturbation, 
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Figure 1: The safety factor defined inside a magnetic island, plotted versus the radial position x, 
and normalized to q(0) = l/(\^Mn is i) = 7/ (\s\ y/1 — e). For typical magnetic islands, 1 — e < 10~ 2 
and we have q(0) > 10. 



Bi s [. The X points of the magnetic island are at {qx — qo, u) = (0, 0) and (0, 2ir) and the O point 
at (0,7r). The flux surfaces of this equilibrium are labeled by ifi = (qr — qo) 2 /2+M(cosu+l). Mis 
a constant determined by the condition • B = (dip/dqT)(dqT/drT)B qT + (dip/du)\ Vu|£? u = 0. 
We obtain M = (qo\s\ /ni S i)(B is i / B po ^ ) , where B po i t0 is the poloidal magnetic field calculated at 
the rational surface. 

2.2 Coordinates inside the magnetic island 

Here, we model the plasma inside the magnetic island as a straight flux tube, where £ is the 
translational symmetry coordinate, and ip labels the nested flux surfaces. The magnetic axis and 
island O point are at ip = 0, while the separatrix is labeled by ip = ip sx = 2M. We define a 
complete set of cylinderlike coordinates (p,9,C), with radial- like and anglelike coordinates given 



With these definitions, the magnetic axis is at p = and the separatrix radius is p = p sx , which 
corresponds to the magnetic island half-width Wi s i, given by the Rutherford formula [21J: 



with e defined below in this paragraph. The angle 8 is defined in the domain (0, tt/2) and extended 
to (0, 7r) by reflection symmetry with respect to 9 = vr/2, with values 0,7r at the rational surface 
q = qo. Further extension to (0, 2tt) is obtained by reflection symmetry for p o —p. We also 
note that the flux surface's cross section in the (p, 9) plane and in the proximity of the O point 
is an ellipse, and the value of its eccentricity is e = 1 — Mn| s >7 2 / ' s 2 . Typical magnetic islands in 
tokamak experiments have eccentricities close to e ~ 1. The case of zero eccentricity corresponds 
to a value of magnetic island width, which is never encountered in tokamak plasmas because a 
saturation or a disruption occurs long before the island grows up to these sizes. Nevertheless, 
in the following sections we calculate the SAW continuous spectrum for the case e ~ too. In 
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fact, this case shows clearly the analogy with the well-known results of continuous spectrum in a 
tokamak, whose cross section has typically small values of eccentricity. The differential operators 
in this coordinate set are provided in Appendix [Aj These are necessary to give an explicit form 
to the equation for the SAW dynamics, that will be introduced in Sec. 13. 11 

2.3 Equilibrium and safety factor 

The equilibrium magnetic field in the region inside the island, can be described in the cylinderlike 
coordinates (p, 9, £) defined in Sec. 12.21 The contravariant physical components are: 

B ph = Q> B e ph = B ph ^P, B< h = B 

Here B® h0 = £oM-Bo/(9o r o7 2 ) is the value of B® h at p = 1, 9 = 7r/2, where all parameters 
refer to the original axisymmetric equilibrium. The function a and the definition of contravariant 
physical components are provided in Appendix |A"1 

We define the safety factor q inside the magnetic-island flux tube as the average of Q = 
pB /(Z B 9 ) over 9, with B e = B°Fy/T=li p); i.e., 



=K(x) = l^J^K(x) (2) 



7r \s\y/l - e ' vr y/Mn is 



where we have introduced the complete elliptic integral of the first kind K(x) = f^ 2 dO/F, with 
F = VT — x 2 cos 2 9 and x = p/p sx - A similar definition for the safety factor was given in Ref. |19j . 
but our result is a factor q$ smaller. The basic difference of our derivation from that of Ref. [19], is 
that we take into account that the length of a magnetic island flux tube is Zq ~ l-nq^R®, whereas 
in Ref. [19] a flux tube with length Zq ~ 2itRq is considered, which is the major circumference of 
the tokamak (see also Ref. [18]). 

The safety factor inside a magnetic island, normalized to the minimum value q{x = 0), is 
shown in Fig. [TJ We see that q(x) has a singular behavior near the separatrix x = 1, since B 8 
vanishes at p = p sx , 9 = 0, ir. This feature is due to the X points and is analogous to safety factor 
profiles in diverted tokamak plasmas. 

3 Shear Alfven wave continuous spectrum 
3.1 Model equations 

Here, we study the shear Alfven wave dynamics in a general equilibrium, adopting the ideal MHD 
model. In next sections we introduce the geometry defined with the coordinates (p, 9, £), and give 
the explicit solution of the problem. The starting linear equations of our model are: 

^ = -VP + - Xj B (3) 

at c L 

E = -— X B (4) 

CL 

dP 

— = - lh PV - v — v ■ VP (5) 

where g is the mass density, v the fluid velocity, P the pressure, J the plasma current, B and E 
the magnetic and electric fields, jh the ratio of the specific heats and cl the speed of light. This set 
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Figure 2: Continuous spectrum of SAW, Q(x), evaluated from the safety factor profile as in 
Eq. \12\ for the case e <C 1 (corresponding to M ~ 1). Modes with n = 10, and m = 9 to 
14 (where m is the poloidal mode number) are considered, and rii s i = 1 is chosen. Here the 
approximation of cylindrical symmetry is adopted, and therefore the continuum branches of the 
modes with different poloidal mode number m intersect without interacting. Compare with Fig. 
and Fig. [H where the continuous spectrum is calculated retaining all geometrical effects due to 
the noncircularity of the flux surface's cross section. 



of equations is solved for the perturbed velocity, the perturbed magnetic field, and the perturbed 
pressure, expressed respectively as 5v = — cx,V<^> X B/B 2 , 5B = V x A and 5p, where <f> is the 
perturbed scalar potential and A is the perturbed vector potential. The set of equations El HI [5] 
for shear Alfven waves can be casted in the following linearized equation for the scalar potential: 
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(j) = 







(6) 



where all quantities but the scalar potential eft, are referred to the equilibrium. All curvature 
effects are retained here: va is a function of space, and the geodesic curvature k s allows the 
pressure to affect the SAW dynamics. The geodesic curvature is defined as the component of the 
curvature n = {B/B) ■ W{B/B) tangent to the magnetic flux surfaces. Equation [6] is valid for 
both global modes and radially localized modes. It can be solved in the coordinate set (p,6,C), 
by substituting the differential operators defined in Appendix lAl 

Here, we want to consider the shear Alfven wave propagation near the resonant flux surfaces 
where the energy is absorbed by continuum damping. Therefore, we focus on the dynamics of 
continuum modes, namely modes that are characterized by radial singular structures. This allows 
us to drop the second and third terms in Eq. The linear equation for radially localized shear 
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Alfven modes in a compressible nonuniform tokamak plasma can be written in the form |14| : 



^-V 2 ± + Z 2 V||ViV||<A - U?BAE - CAP V\4> = (7) 



where uja = va/Zq. We adopt the value of the frequency of the low- frequency SAW continuum 
accumulation point (CAP), delimiting the frequency gap of the beta induced Alfven eigenmode 
(BAE) HHI51, given in Ref. JT6] : 



1 2T i( 7 T e . 
UBAE-CAP = S~\ - + — (8) 



R V rrii V4 T, 

where Tj and T e are the ion and electron temperatures, and rrii is the ion mass. Here, we focus 
on frequencies higher than ujbae-cap and consistently neglect kinetic effects associated with 
wave-particle resonances [16] . The operators Vy and are the gradients calculated along and 
perpendicularly to the equilibrium magnetic field (given by the sum of the axisymmetric field 
plus the island field). 

3.2 Eigenvalue problem 

In section 12.21 we have described the region inside a magnetic- island flux tube with a cylinderlike 
set of coordinates, labeling the magnetic flux surfaces with a radial-like coordinate p. This 
magnetic flux tube is assumed here to be the static equilibrium of the problem. We want to 
put the equation for radially localized shear Alfven modes, Eq. which is written in a general 
vectorial form, in an explicit form, using a cylinderlike set of coordinates for a straight flux tube 
and the differential operators defined in Appendix[A] We substitute the parallel and perpendicular 
differential operators: 

Z V|| = ^Mn isl F— + — (9) 

pI^I = a ^2 ( 10 ) 

and we impose boundary conditions of periodicity in 9 and ( on 2ir circles. Moreover, using 
the symmetry of the equilibrium in the ( coordinate, we can write the general solution as <ft = 
4>(x, 9) exp (— inC), where n is the mode number in the ("-direction. 

We also note that the problem we are considering is greatly simplified if we change coordinates 
from (p,9,C) to a field aligned set of coordinates (p,9,^), where £ = £ — x( x i^), and x = 
(^fMnisi)^ 1 J® d9'/F(x, 9') is the normalized incomplete elliptic integral of the first kind. In fact, 
we have B ■ V£ = and the parallel derivative takes the simple form ZqV\\ = \fMrii s iF d/d9. 
Finally, Eq. [7J is written in the form of an eigenvalue problem: 



9? F d d 
+ — -— aF- 



iMn 2 , a 86 89 

isl 



(11) 



where Q 2 = (uj 2 — ^bae-cap) / u a ^ s ^ ne eigenvalue, and <fi" = (8 2 c/)/d 2 x) exp (— inx) is the 
eigenf unction. The boundary condition <j)"(9 = 0) = <p"(9 = 2ir) now reads (j)"{9 = 0) = <p"(9 = 
2ir) exp (2mnq), with the safety factor q defined in Eq. [2J 
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Figure 3: Continuous spectrum of SAW, Q(x), obtained by solving numerically Eq. [TT], and 
plotted versus the radial position inside the island, x. The case of small eccentricity case, e<l 
(corresponding toM~ 1), is considered. Typical values of the equilibrium parameters have been 
chosen and rii s i = 1. The O point is at x = and the separatrix at x = 1. Modes with n = 10 
are considered. The MiAE gap is found at frequencies Vl ~ \fMrii s i. 



Before solving Eq.[TT]in its exact form, we note here that we can have an approximate solution 
by assuming cylindrical symmetry, and using the safety factor as calculated in Eq. [2j In this case, 
the solution of Eq. [TT] can be written in the form: 

Q = {nq — m)/q (12) 

where we have used d/dO = im, with m being the poloidal mode number. When cylindrical 
symmetry is assumed, the continuum branches of modes with different m are not coupled, and 
therefore intersect without interacting (see Fig. [2]). In the case of small eccentricity e <C 1, and 
close to the O point x ~ 0, we have a = F = 1 and the formula given in Eq. [12] is a good 
approximation of the solution of Eq. [TT] ft = y/Mn,i S i(nq — m). 

In next section, we provide the exact solution of Eq. [TT| namely retaining all geometrical 
effects due to the noncircularity of the flux surface's cross section. 

3.3 Solution 

In previous section, the equation for SAW continuous spectrum has been written in the form of an 
eigenvalue problem in the coordinates (p, 9, £), Eq. [TT] This is an ordinary differential equation in 
one variable, 9, where the radial position x is treated as a parameter. It can be solved numerically, 
with a shooting method code for each position < x < 1, giving the continuous spectrum Q 2 (x) 
as result. 

The case of small eccentricity, corresponding to M = 1 (and e <C 1), is considered firstly. 
Typical values for the equilibrium parameters, qo = 2, s = 1, £o = 0.1, = 1 are chosen. 
Continuum modes with mode number n = 10 are considered (see Fig. [3]). In this case of small 
eccentricity, we see that near the O point, x ~ 0, cylindrical symmetry is a good approximation 
and the result is well described by the formula given in Eq. [T2] as shown in Fig. [2] On the 
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Figure 4: Continuous spectrum of SAW, fi(cc), calculated as in Fig. [3J for n = 1 and a typical 
size magnetic island, M = 1CP 2 (corresponding to e ~ 0.99). In this case, the structure of the 
continuous spectrum is the same as for M = 1, but the MiAE gap is much wider. This implies 
that a wide range of frequencies exists, where modes contained within the magnetic island are 
not affected by continuum damping. 



other hand, for x ~ 1, the flux surfaces have a noncircular cross section and therefore mode 
coupling occurs between modes with different m numbers. This mode coupling is responsible for 
the presence of a gap in the spectrum analogous to the EAE gap in tokamaks [20J, dubbed here 
as magnetic-island induced AE (MiAE) gap. The closer a flux surface is to the separatrix, the 
wider is the MiAE gap. 

The calculation is also repeated for a magnetic island with a typical size, with M = 1CP 2 , 
using the same equilibrium parameters. This corresponds to a magnetic island with Wi s i ~ 
\/l — efQ = 0.1 ro, which is the order of magnitude of a saturated magnetic island in tokamaks. 
For this case, we considered modes with n = 1. We find that the continuous spectrum structure, 
shown in Fig. is the same as for M = 1, but in this case we obtain a much wider MiAE gap. 
This is due to the fact that flux surfaces of a typical magnetic island have a high eccentricity, and 
therefore the mode coupling between modes with different m numbers is stronger. The MiAE-gap 
central frequency is proportional to the magnetic island half-width and near the O point we have: 

n nr7 qosn is iWisi 

SlMiAE = VMn is i = — (13) 

2 r 

This is valid for any value of any value of magnetic island half width, corresponding to any value 
of eccentricity. 

For typical tokamak plasma parameters, u bae _ cap /lu a ~ 0Qoi w hh (3 denoting the ratio 
between plasma and magnetic pressures. Therefore, finite compression BAE frequency gap and 
magnetic island induced MiAE gap are of comparable size when M ~ ftq^, i.e. for a typical 
saturated magnetic-island size in tokamaks. We have considered a straight flux tube model, 
neglecting the effect of curvature coupling among different m numbers. Therefore, a gap analogous 
to the toroidicity induced Alfven eigenmode gap [HJ [22] is not found here. Nevertheless, such a 
toroidal MiAE gap is expected to be present inside magnetic islands, at frequencies lower than 
the frequencies of the ellipticity MiAE gap. Its gap width can be estimated by substituting va 
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with va/ (1 + £q cos 9t) in Eq. [71 This gives a multiplication factor (1 + 2eq cos 9t) to the first and 
third terms of Eq. and consequently to the first term of Eq. [TTJ Therefore, we expect to have 
a toroidicity induced gap within a magnetic island, with a width of the order of AO 2 ~ Mnf sl eo, 
analogously to the theory of toroidicity induced AE. The magnetic island also has a modulation 
in the tokamak toroidal direction, given by the number n^. This helicity of the flux tube is 
neglected here and is expected to weakly couple modes of the SAW continuous spectrum with 
different n, analogously to the case of stellarators. These effects are not expected to qualitatively 
modify the MiAE gap and will be investigated in a different work as further extension of the 
present theory. On the other hand, the BAE gap is considered in our model in the definition of 
il 2 , resulting in an upward shift in the continuous spectrum of the frequencies uj 2 /uj\. 

4 Conclusions and discussion 

One of the main damping mechanisms of SAW instabilities in nonuniform plasmas is continuum 
damping, which occurs near the magnetic flux surfaces where the frequency of the instability 
matches the SAW continuum frequency. Therefore, understanding the structure of the SAW 
continuous spectrum is one crucial step in studying the stability properties of a tokamak plasma. 
In toroidal devices, a frequency gap in the continuous spectrum represents a range of frequency 
where discrete Alfven eigenmodes (AE) can grow unstable, practically unaffected by continuum 
damping. 

The radial structure of the continuous spectrum of SAW is calculated in this paper in an 
equilibrium inside a magnetic-island flux tube. Since the typical island frequency and growth 
rate are much lower than the SAW oscillation frequency, we have assumed that the equilibrium 
magnetic field is given by the usual tokamak axisymmetric field plus a quasi-static helical distor- 
tion due to the magnetic island. The flux tube is considered to be straight, in order to focus on 
the effects of the noncircularity of the cross section. Curvature effects are retained only in the 
formation of the BAE gap in the low frequency part of the continuous spectrum. A linear MHD 
model is adopted for finite-beta tokamak plasmas. A generalized safety factor is defined inside 
the magnetic-island flux tube, giving information on the rational flux surfaces where different 
modes couple to create gaps in the SAW continuous spectrum. The linearized SAW equation is 
solved with a shooting method code inside a finite-size magnetic island and the result compared 
with the SAW continuous spectrum calculated in tokamak equilibria. 

We have found that there exists a SAW continuous spectrum within a magnetic island, similar 
to that calculated in tokamak equilibria [TTJ. For the case of small eccentricity of the flux sur- 
faces, modes with different poloidal numbers are not coupled and therefore, the branches of the 
continuous spectrum intersect and there is no formation of gaps. This is shown to occur inside 
a magnetic island with very large (not typical) width, near the O point, where the flux surfaces 
are nearly circular. On the other hand, a typical-size magnetic island is shown to have wide gaps 
in the continuous spectrum, due to the strong eccentricity of the flux surfaces, labeled here as 
MiAE gap. This is analogous to the formation of EAE gaps in tokamaks, but is found here inside 
the magnetic island equilibrium. Note that MiAE can exist as bound states within the island, 
essentially free of continuum damping, provided that plasma equilibrium effects and free energy 
sources can drive and bind them locally [1]. 

The equilibrium inside a magnetic island is shown here to be analogous to a tokamak equi- 
librium, and, consequently, we emphasize the analogies between AE within the magnetic islands 
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(MiAE) and the well-known AE in tokamaks. The MiAE peculiarity resides in the corresponding 
SAW continuous spectrum dependence on the magnetic-island size. Moreover, we expect to find 
all kinds of Alfvenic instabilities inside a magnetic island, and not only ellipticity-MiAE. For 
instance, toroidicity-MiAE could grow inside a gap due to the flux tube toroidicity effects, which 
are not studied here and will be treated in a subsequent paper. 

We point out the important implications of this result in understanding the dynamics and 
stability properties of a magnetic island. In fact, the presence of MiAE inside a magnetic island 
could modify the equilibrium profiles and nonlinearly affect the magnetic island growth. On 
the other hand, MiAE could nonlinearly affect energetic particles redistribution in the proximity 
of the magnetic island rational surface, caused in part by the radial magnetic field due to the 
magnetic island itself [23] . 

The frequency of the ellipticity MiAE gap has been shown to be proportional to the magnetic- 
island width [see Eq. [J3] . This suggests the possibility of using the MiAE frequency scalings as a 
novel magnetic-island diagnostics in a fashion similar to other commonly used Alfven spectroscopy 
techniques |244 I2"5] 126] . Unlike in the BAE case, the radial MiAE localization at the center 
of the island makes them difficult to detect by external measurements, and makes the use of 
internal fluctuation diagnostics, such as electron cyclotron emission and soft x rays, necessary. 
By detecting MiAE inside the island and measuring their frequency, one has indirect information 
on the magnetic-island size. 
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A Differential operators 



Here, we calculate the differential operators in the set of coordinates inside the magnetic island 
(p, 9,Q- The domains of these coordinates are < p < p sx = Wi s i, < 9 < 2ir and < ( < 2ir, 
where ( is the coordinate of translational symmetry, and periodicity in 9 and £ is assumed for 
the perturbations. The gradients of these coordinates are: 



V/9 = sin 9 q — cos 9\J\ — eF u = \fa p 

V(9 = (cos#g + sin<Vl - eFu)/p = VbO/p 

vc = t/z 

where we use the notation V = V /V for a unitary length vector. Moreover we have the relation 
Vp • V# = c/p. Here F = a/1 — x 2 cos 2 9, where x = p/p sx , and: 

a = sin 2 9 + cos 2 9(1 - e)F 2 
b = cos 2 9 + sm 2 9(l-e)F 2 
c = cosflsinfl(l - (1 - e)F 2 ) 

A covariant basis (g P ,g9,g() can also be defined as orthogonal to the contravariant basis 
(g p = V/9, g 6 = V#,g^ = VC), namely satisfying: g % gj = 5*. With this notation, any vector can 
be written as V = J2i V % gi = Yli Vph9i> where we call V^ h the physical components of the vector, 
Vp h = V l \gi\. Finally, the differential operators for the cylinderlike set of coordinates (r, 9, £) are: 

v/ = ^ + (14) 

op p 39 Zq o( 

v v ~ z d( + p Up I Vb ) + de\^)) (15) 

V/ " pdp\FVdp + pd9)) + pd9\FVdp + pd9)) + Z 2 d( 2 ( ' 

_ Vbdv P i i / c dv; h ^v- b dv° h] 

(VXVjp " p 89 + FVT^\Z d( Z d( ) {U) 

m v V x r dv P l + 1 ( ^ dv & c dv h \ 

(VXV)9 = + -F-JT=e\—^K - Z- lK) (18) 



1 jC pVdpV fV& ^ > d0\ F ^ F^EJJ 



(19) 



It is straightforward to see that, in the case e ~ and in proximity of the O point (F = a = b = 1, 
c = 0), these differential operators reduce to the known differential operators in cylinder geometry. 
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